
I Office J 



rci/bb ZUU4 / U U 2 9 U << 




INVESTOR IN PEOPLE 

The Patent Office 
Concept House 
Cardiff Road 
Newport 
South Wales 
NP1Q 



m§ 

recd 2 3 JUL 200*1 



WIPO 



PCT 



I, the undersigned, being an officer duly authorised in accordance with Section 74(1) and (4) 
of the Deregulation & Contracting Out Act 1994, to sign and issue certificates on behalf of the 
Comptroller-General, hereby certify that annexed hereto is a true copy of the documents as 
originally filed in connection with the patent application identified therein. 



In accordance with the Patents (Companies Re-registration) Rules 1982, if a company named 

in this certificate and any accompanying documents has re-registered under the Companies Act y=J 

1980 with the same name as that with which it was registered immediately before re- 

registration save for the substitution as, or inclusion as, the last part of the name of the words BJLE 

"public limited company" or their equivalents in Welsh, references to the name of the company 

in this certificate and any accompanying documents shall be treated as references to the name 

with which it is so re-registered. 

In accordance with the rules, the words "public limited company" may be replaced by p. I.e. , ^ ^ 

pic, P.L,C. or PLC. ^ 

Re-registration under the Companies Act does not constitute a new legal entity but merely 
subjects the company to certain additional company law rules. 




patents Form 1/77 

k.ctl977 



^Patents 



Request for grant of a patent 

(See the notes on the back of this form. You can also get 
an expUnatory leaflet from tlte Patent Office to help 
you fill m this form) JJ^wneip 



2. Patent application number 
(The Patent Office w ill fiU m this part) 

3. FuU name, address and postcode of the or of 
each applicant (underline all surnames) 



It) 



Patents ADP number (ifyou know 

If the applicant is a corporate body, give the 
country/state of its incorporation 



4. Tide of the invention 




IP/P7295 



/ / / 



The Patent Office 

Cardiff Road 
Newport 
Gwent NP9 1RH 



0316402.7 



P01/7700 0.00-0316402.? 



Registered Office 85 Buckingham Gate 
London SW1E 6PD 
United Kingdom 



GB 



5. Name of your agent (if you have one) 

"Address for service" in the United Kingdom 
to which all correspondence should be sent 
(including the postcode) 



Stents ADP number (ifyou kno w it) 

11 you are declaring priority from one or more 



Bowdery, Anthony Oliver 

QINETIQ LIMITED 
IP Formalities 
A4 Bldg 

Cody Technology Park 

IvelyRoad 

Farnborough 

Hants GU14 0LX United Kingdom 



earlier patent applications, give the country 
and the date of filing of the or of each of these 
earlier applications and (ifyou know it) the or 
each app lication number 

If this application is divided or otherwise 
derived from an earlier UK application, 
give the number and the filing date of 
the earlier application 



Country 



rionty application number 
(ifyou know it) 



Number or earlier application 



Date of filing 

(day /month /year) 



Date of filing 
(day /month /year) 



Is a statement of inventorship and of right 
to grant of a patent required in support of 
this request? (Answer •Yes' if: 

t Z.T ifca " f ' UWied ln P art ^>u>tan inventor, or 
o) mere is an inventor who is not named as an 
applicant, or 

c) ™y named applicant is a corporate body 
See note (d)) y 



YES (b) 



itents Form 1/77 



Patents Form 1/77 



Enter the number of sheets for any of the 
flowing items you are filing with this form. 
x/o not count copies of the same document 

Continuation sheets of this form 0 

Description 13 

Claims 03 

Abstract 01 

Drawing^ 



;t 01 v 




10. If you are also filing any of the following, 
state how many against each item. 

Priority documents 0 

Translations of priority documents 0 

Statement of inventorship and right 01 
to grant of a patent (Patents Form 7/77) 

Request for preliminary examination 01 
and search (Patents Form 9/77) 

Request for substantive examination 0 
(Patents Form 10/77) 

Any other documents 0 
(please specify) 



11. I /We request the grant of a patent on the basis of this application. 

Signature J^T"^^^^^ — ~r 



Date 11 July 2003 
Mr A W S Williams 



12. Name and daytime telephone number of 

person to contact in the United Kingdom Mrs Linda Bruckshaw 01252 392722 

Warning __ 

After an application for a patent has been filed, the Comptroller of the Patent Office will consider whether publication 
or communication of the invention should be prohibited or restricted under Section 22 of the Patents Act 1977. You 
will be informed if it is necessary to prohibit or restrict your invention in this way. Furthermore, if you live in the 
United Kingdom, Section 23 of the Patents Act 1977 stops you from applying for a patent abroad without first getting 
written permission from the Patent Office unless an application has been filed at least 6 weeks beforehand in the 
United Kingdom for a patent of the same invention and either no direction prohibiting publication or communication 
lias been given, or any such direction has been revoked. 

Notes 

a) If you need help to fill in this form or have any questions, please contact the Patent Office on 0645 500505. 

b) Write your answers in capital letters using black ink or you may type them. 

c) If there is not enough space for all the relevant details on any part of this form, please continue on a separate 
sheet of paper and write "see continuation sheet" in the relevant part(s). Any continuation sheet should be attached 
to this form 

d) If you liave attached 'Yes* Patents Form 7/77 will need to be filed. 

e) Once you have filled in the form you must remember to sign and date it 

f) For details of the fee and ways to pay please contact the Patent Office. 



Patents Form 1/77 



1 



Direction Finding 



10 



This invention relates to direction finding, and more particularly to direction finding using 
radio techniques. 

Direction finding by radio is known. In a direction finding system having a number of 
separate antennas, radio wavefronts reach the antennas with different delays. Assuming 
the wavefronts have narrow bandwidth, these delays give rise solely to relative phase shifts. 
For determining the directions or bearings of M-1 radio frequency (RF) emitters (M is a 
positive integer), a typical direction finding system employs M antennas, a respective 
receiver for each antenna, and processing circuitry for each receiver. The processing 
circuitry implements a discrete Fourier transform to divide its respective signal into 
frequency bins or channels. Signals are then combined in pairs, and for each channel a 
respective covariance matrix is constructed from which emitter bearings are estimated." This 
approach suffers from the problem of requiring as many receivers as there are antennas, 
receivers being expensive and bulky. It is unsuitable for example for man-portable 
15 equipment or for mounting in small aircraft. 

To get round the multiple receiver problem, it has been proposed to use a single port 
receiver, spatial spectrum estimation technology, and a weight perturbation algorithm to 
obtain the covariance matrix. See Zhao Yimin, "A Single Port Receiver Spatial Spectrum 
Estimate DF System", 0-7803-3216-4/96 IEEE 1996. This is however a relatively complex 
20 approach to the problem. 

It has also been proposed to combine RF antenna signals in beamformers and use a single 
receiver connected to successive beamformers via a multipole switch, one pole per 
beamformer. See CMS See, "High Resolution DF with a Single Channel Receiver", 0- 
7803-7011-2/01 IEEE 2001. A single spatial covariance matrix is formed, and the bearings 
or directions of up to M-1 emitters can be estimated, where M is the number of antennas. 
This requires M 2 beamformers, another source of expense and bulk. 

Another solution is adopted in the Rohde & Shwartz DDF 195 instrument, which combines 
pairs of RF antenna signals with each of four relative phase shifts inserted between them in 
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succession by means of switches This requires only a single receiver, and uses one 
antenna as a reference antenna, combining its output with that of other antennas in turn 
with multiple switched phase shifts. However, the method is for estimating the bearing of a 
single emitter only. A patents search has indicated that the following patent documents 
EP455102, DE4014407, DE3636630, DE1 9529271 and DE2723746 are related to direction 
finding. 

US patent application no. US 2002/0190902 A1 describes sampling RF signals from M 
antennas at a sampling rate equal to or greater than the signal bandwidth multiplied by 2M. 
This is followed by Fourier transformation of resulting signal samples to provide spectra. 
Direction finding is then based on line configuration in the spectra and associated phase 
and amplitude data. It employs a multiplicity of directional antennas to provide 
beamforming. 

It is an object of the present invention to provide an alternative form of direction finding. 

The present invention provides a direction finding system incorporating a plurality of 
antennas characterised in that the system also includes: 

a) means for determining individual antenna signal strengths; 

b) combining means for determining combined antenna signal strengths by forming 
combinations of first and second antenna signals derived from different antennas, 
wherein the second antenna signals are in two sets with signals in one set having a 
non-zero phase difference relative to signals the other set; and 

c) means for deriving covariance matrix elements from antenna signal strengths and for 
determining at least one emitter bearing therefrom. 

The invention provides the advantage of requiring only a single receiver when successive 
signal strengths are determined in successive steps, which reduces cost and bulk and still 
provides a viable direction finding technique. 

The relative phase difference may be in the range 30 to 120 degrees, preferably 
substantially 90 degrees. 
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The combining means may be arranged to -combine antenna signals with equal gain 
magnitude and with equal and unequal phase. It may incorporate phase shifting means 
switchable into and out of an antenna signal path, and an adder having two inputs both 
switchably connected to individual signal paths extending to respective antennas. 

The means for determining individual antenna signal strengths may comprise a first 
multipole switch having input poles connected to receive signals from respective antennas; 
the combining means may incorporate a second multipole switch having input poles 
connected to receive signals from respective antennas and a third multipole switch for 
switching phase shifting means into and out of an antenna signal path extending via the 
second multipole switch; and the combining means may also incorporate adding means for 
combining signals, the adding means being arranged to add an antenna signal in a first 
signal path extending via the first multipole switch to another, antenna signal in a second 
signal path extending via the second and third multipole switches. 

In another aspect, the present invention provides a method of direction finding using a 
15 plurality of antennas characterised in that the method incorporates: 

a) determining individual antenna signal strengths; ' 

b) determining combined antenna signal strengths by forming combinations of first and 
second antenna signals derived from different antennas, wherein the second antenna 
signals are in two sets with signals in one set having a non-zero phase difference 

20 relative to signals the other set; and 

c) deriving covariance matrix elements from antenna signal strengths and determining at 
least one emitter bearing therefrom. 

The relative phase difference may be in the range 30 to 120 degrees, preferably 
substantially 90 degrees, and successive signal strengths may be determined in successive 
25 steps. 

The step of forming combined antenna signal strengths combines antenna signals with 
equal gain magnitude and with equal and unequal phase. It may include switching phase 




shifting means into and out of an antenna signal path, and adding signals In signal paths 
extending switohably to different antennas. 

The step o, determining IndMdual antenna signal strengths may comprise swHching signais 

a) switching signals from antennas via a first path for combining; 

b) switching signals from antennae via a switch selectable second path or a third path for 
comb,n,ng, the third pafh being arranged to phase shift antenna signals therein relative 
to antenna signals in the second path; and 

10 c) adding a first path antenna signa to second and third pafh antenna signals Individually. 

In order the, the invention might be more fully understood, embodiments of the Invention will 

inlteh ^ ^ ^ °' eXamP ' e ° nly ' re,erenCS t0 ,h6 -~™ 9 twinge, 
15 Figure 1 is a schematic block diagram of a direcSon finding system of the invention; 

Figure 2 is a graph illustrating direction finding results obtained using the invention; 

Figure 3. is a graph Illustrating direction finding results obtained using a single receiver 
pnor art system comparable with the invention; and 

Figure 4 is a graph lliustratlng direcion finding results obtained using a multiple receiver 
zu prior art system; 

Referring to Rgure 1, a direction finding system of the Invention is indicated generally by 
10 As „,us.ra,ed for .he purposes o, this example, the system to incorporates four 
antennae 12 eaoh * whloh ls . g 

antennas may be used as are required to detec a desired number of emitters, i.e M 
25 antennae for M-1 emtoers. Signais pass from me antennas 12 vta respective buffer 
amplifiers ,14 to firs, and second multipole switches SW1 end SW2, the amplifiers 12 being 

ziTzz'r*" ,nput b ' °' d ^ e - b ° ,h ^ 

SW and SW2 have respective movable confacte 01 , C2 which allow any of .he associated 
.nput poles b to d in each case to be connected to respective oufpu. poles f. 



The outpu po,e f o> the flrs, switch SW1 is connected «o a firs, phase sh«er P1 inducing 
a phase 8h « of and ths output pQ|e , of ^ ^ ^ ^ ^ ^ 8 

pin ell r P2 ^ P3 ' n,r0dUC ' n9 PhaM **• * * and * «~ely. in the 
present exampie, where an ideal situation is envisaged, ail three phase sheers P1 ,o P3 

have a ga,n of unity: the flrs. and second phase shifts and are zero (and flte 

associated phase shifters P1 and P2 couid he removed and replaced by connection; 2 

■ 1 h T ^ " 90 de9 ' eeS - H ° WeVer ' n °"- Wsal »aees can be 

ZT T y Ca " bra " 0n <" ^ deS ° ribed ^ Sa,i ^ P^-noe he 
een demonstrated wflh errors up to plus or minus 60 degrees: in omer words the 
10 d fterence * deduced by *. third phas6 shKer p3 may fae ^ ^ ^ ^ . 

± 20 degrees, with simultaneous gain discrepancy of up ,o 3 dB compared to unity. Larger 
sorep ,„ progressive deten . orafcn of performance ^ 

I PrOVld ' nS * ,0 ^ *» ™" 30 .o ,20 degrees different to 

15 provl: TT*" " * a ' S0 aC06P,ab ' e ,0f PhaSS Sm a " d 9 ai " '° «» * <~y, 
15 p ro „ded the phase shift is reasonably Cose to 90 degrees and phase shifter ga ns are 

"eToy ,0 1 °' ^ ,ha ' ^ ~° ™ ^ - P^-rmined funlnl 

^u!s W T fr0m PhaSS P1 PaSSSS ,0 3 « *" 16a - - adder 18 

20 npu, pol s p and q respectively o, a third multipole switch SW3, which has a third input pole 
~ 8 '* * ird 3 m ° Vabte — C3 «*■* any 

* nn r t r p ; les * q and r ,o be ,o « -*« p ° ie * * - - 

connected to a second input 1 6b of the adder 1 8. 

25 perfect I'd IT"*" draW ' n8 "** a " °° m ™ ™ assumed to be 

may need to be fnmmed or impedance matched to counteract unwanted effecte by inserting 
addtflona, crcuitry. This is we,, known ,n fhe art o, electronics and win no, be described 



30 



Isla *** 20 '° 3 aaHee-ocnneofed ohain of eiemente 

Zc mer Is TT ' " * C ™ (AD °> 24 ' a *«* Fourier 

transformer 26 and a dtgrtal signal processor (DSP) 28. The digital elements 24 to 28 may 



touted by equivalent analogue processing , de s.red. The transfer 26 provides a 
derate Founer transform (DFT, of .ength N points with windowing. This is a well-known 

The : 6 Si9na ' Pr ° CeSSin9 "' "* En9teWOOd "SA. « 

~ 'Theth ^ ' aSt F ° Urier ,ranS, ° m <FFT) ' b * ,hls te * 
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I, L h ' 9 SyS,6m 10 ° Pera,eS 38 ,O " 0WS - By appr ° pria,e °< *e positions 

o, he movable sw„ch contacts C, ,o 03, any antenna signa, may be added ,o any other 
antenna srgna,, or seieoted aiona by adding i, ,o a zero signa, a. grounded inpu, r « 
swrtch SW3. Assuming ma, the three phase shifters Pi to P3 have aqua, gain magn^es 
antenna srgnais a, the adder inputs 16a and ,6b have a reiativa phase sh« be twee^em 
oqua, to or W aocording to whether the third switch movabie oontact 03 Is 
connected to fts fire, or seoond input pole p or q respectively. In ma ideal case, ♦,.«,.„ 

bale = T 9 " 8 ' *"* *" re,at ' Ve SW " ° °' 90 <'gno„ng sign, 

between adder ,npu, signals according to whether the third switch movable 

connected to its first or seoond input pole p or q. 

olTaT ^ ^ V**" " fi ' ter (n0t Sh ° Wn> '° ~ e »™<*« ^na,s 

ou s,de a frequency band in which the antenna signals appear, ft processes the adder 

.put s.gna, and converts ft to a compiax output i„ a base band convenient for sampllg 

and analogue to digital conversion a, 24. The sampHng rafe o, the ADC 24 J 2 

uff,c,en,,y n ,gh to avoid aiiasing, i.e. i, must exceed the bandwidth o, the receiver Zl d 

bandpass After, if, for example, ma fronf-and bandwidth is 25 MHz, and a 512-poin Dm 

mpiemented by ,he transformer 26, ,he ADC sampling rate would be a. least 27.5MHz and 

ft wouid take approximately 20 pa to acquire 512 samptes for the transformer 26. The digfta, 

MreJr ^ * °° nVerted ^ ^ ,ranS,0rmer 26 * a *~ — - • ~ 

1 nlde ir"^ M " Channe ' S °' *» *' h SaCh ^ » — — 
DFT vlth The | . ,reqU6n0y ° hanne,S WOU,d be approximately 100 kHz wide for a 512 point 

bal k 7 27 ' 5MHZ - A IS imP ' emented Wi * *> -dues ieakage 

112,°, 0P,i ° nS Ch ° iCe ° f Wind ° W * *dow 

equivalent to no window. 



7 



Optionally, selection of settings of movable contacts C1 to C3 may be made in such a way 
that not all antennas contribute. This gives faster processing, but possibly less accuracy, 
and it reduces the maximum number of emitters that can be detected. 

As has been mentioned, the various possible settings of ihe switches SW1, SW2 and SW3 
5 allow the receiver 22 to input either the signal from any individual antenna 12, or a sum of 
relatively phase shifted signals from any pair of antennas 12. The switches SW1 to SW3 
may be operated to give a random or pseudo-random selection of antenna signals to avoid 
possible deleterious effects with particular signals. The output of the receiver 22 is sampled 
■ by the ADC 24 and processed by the transformer 26, which computes the N-point 

10 windowed DFT of a block of N consecutive samples from the ADC. This is a well known 
procedure and will not be described in detail: The output of the transformer 26 comprises N 
frequency domain samples, i.e. magnitudes of the contents of the N frequency bins. For the 
general or nth frequency bin, the frequency domain sample is denoted by S n , where n is a 
frequency domain index in the range 0 to -N-1). sf is defined as the discrete Fourier 

15 transformer output with frequency index n if the kth antenna is connected through the first 
swtch SW1, the third switch SW3 is connected to OV at r, and the gain of the entire path 
from the kth antenna to the receiver input is unity. 

The process of direction finding consists of estimating the bearing or angle of incidence or 
angle of arrival of one or more signals received by the antennas 12. It is carried out for one 
20 or more frequency bins by the digital signal processor (DSP) 28, which processes frequency 
domain samples S n . 

A complex gain constant G 1k is now defined to represent both gain and phase shift applied 
to an antenna signal in a path from the kth antenna 12 (k = 1 to M). through the first switch 
SW1 , the first phase shifter P1 and the adder 20 to the receiver 22. If the buffer amplifiers 
14 have equal gains, and input poles b to d of the first switch SW1 have like properties, 
then G 1k has the same value for all antennas 12 and is written G,. If the gains of the buffer 
amplifiers 14 are not equal, or the first switch SW1 properties vary between inputs, the 
value of G 1k will be different for each antenna 12. To simply analysis, these gains are 
assumed to be equal for all antennas 12: If they are not in fact equal, the term G, should be 
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replaced by G nk in the analysis below, which is then carried out for each antenna 12 
separately. 

Similarly, complex gain constants and G^ are now defined to represent both gain and 
phase shift applied to an antenna signal in paths to the receiver 22 from the kth antenna 12 
via (inter alia) the second and third phase shifters P2 and P3 respectively. Again, G^ and 
Ga* are assumed to be the same for all antennas. 12, .and are written G^ and G 2B 
respectively. As before, if this assumption is untrue, the gains should be replaced by G^ 
and Gsbk and the analysis carried out for each antenna 12. 

The gain of the combination of the receiver 22, ADC 24 and transformer 26 do not affect the 
direction finding result, so to simplify analysis without affecting results this gain is assumed 
to be unity. 

The gains G n , G^ and G^ are first measured at even/ frequency of interest, i.e. at the 
centre frequencies of those frequency bins defined by the DFT operation in the transformer 
26 which are associated with emitters. With the first switch SW1 connected to the buffer 
amplifier 14 of the kth antenna and the third switch SW3 connected to OV at r, the 
transformer output with frequency index n is the product Q, S£ . This transformer output has 
a mean squared value or power P kk associated with the kth antenna and given by: 

_P kk =E{|G 1 5 fe "| 2 } = |G l | 2 E{|^| 2 } (1) 

where E{..} is the mathematical expectation operator and |...| (as in e.g. |S« |) represents a 
modulus. The DSP 28 computes a measurement (or "estimate") of P kk which equals either 
a single value of S n or an average or weighted average of several values of S n , these 
values being obtained from the discrete Fourier transforms of respective blocks of data 
collected with switch contacts C1 to C3 set to appropriate positions. These blocks of data 
may be overlapping or non-overlapping. For clarity in the following explanation, P kk will 
25 denote the measurement of P kk obtained as described above. 
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The first switch SW1 is connected to the buffer amplifier 14 of each antenna 12 in turn, i.e. 
the antenna index k goes from 1 to M where M is the number of antennas, and the third 
switch SW3 remains connected to 0V at r. The transformer output power P* is measured in 
each case. E{|S£ | 2 } is then computed for each antenna 12 using the previously measured 
value of Gi as: 

E{|^| 2 }=P,*/|G 1 | 2 ' (2) . 

The first switch SW1 is now connected to the kth antenna, the second switch SW2 is 
connected to the mth antenna, and the second switch signal path with gain G^ is selected 
by connecting the third switch movable contact C3 to its first input p. The power P kmA 
associated with gain G^ at the transformer output with frequency bin index n is then 
measured, and it is given by: 

PkmA = £{10,^ + G 2A S£| 2 } (3) 
i.e. P,™ = IG^Efl S>{ | 2 } + | G2A | 2 E{| S£ | 2 } + 2Re{ G^Ef S^S" * }} (4) 

where Re{..} represents "real part of" and the asterisk "*" a complex conjugate. P,^ is 
measured for all pairs of different antennas 12, i.e. for antenna index k = 1 to M-1 paired 
with m = k+1 to M respectively. This implies an antenna is not paired with itself. The value 
of P^ is measured by the DSP 28 in the same way as it measures P,*, as previously 
described, and the same convention is adopted that the' notation P^ is used in what 
follows to refer to the measurement. 

20 Using the measured values of G 1f G^ and the values of E{|,S£| 2 } and E{|S£| 2 } computed 
above, a quantity x is now computed for each antenna pairing from: 

x = d/2)(P kmA - IG^Efl S£ | 2 } - | G2A | 2 E{| S» | 2 }) . (5) 
i.e. x = Re{G 1 G* 2A E{S£S£*}} . (6) 
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In the same way, the procedure associated with Equations (3) to (6) is now repeated, 
except that the second switch signal path with gain G 2 b is now selected by connecting the 
third switch movable contact C3 to its second input q to implement phase shift <t> 3 instead of 
<t> 2 . The power PkmB associated with gain G2B at the transformer output with frequency index 
n is then measured for each antenna pairing, and equivalents of Equations (3) to (6) may be 
generated by replacing index A with index B. This enables a quantity y to be calculated from 
equivalents of Equations (5) and (6): 

y = (1/2)^ - |G n | 2 E{| S£ | 2 } - |G2b| 2 E{| S» ft) (7) 
i.e. y = Re{ G,G\ B E{SJ^S^ *}} (8) 

The known complex value G 1 G* 2A is now written as c + jd, and that of Gn G*^ ase + jf, 
where j is the square root of -1 . The next step is to compute the unknown complex value 

E i s k s m * } written as a + jb. Rewriting Equations (6) and (8) in terms of a to e and j: 

x = Re{(c+jd)(a+jb)} = ca - db (g) 

and y = Re{(e+jf)(a+jb)} = ea - fb (1 0) 

Equations (9) and (10) are two simultaneous equations in two unknowns which are solved 
by standard methods to give the required values a and b, which in turn give E{S^S^ *}. 

For example, for G, = 1 , G2A = 1 and G 2B = j for <j) = 90 degrees, then c = 1,d = 0,e = 0 and 
f =-1, and a = x, b = -y. 

The procedure associated with Equations (3) to (10) is repeated for each chosen pair of 
antenna index values k and m, and the corresponding value of E{Sj^S^ *} is computed in 
each case. These two procedures yield a set of values EflSX) 2 } for k = 1 to M and 
^i s k s m *} for k = 1 to M-1 and m = k+1 to M which are known collectively as "covariance 
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terms". When arranged in a square array with E{| S% | 2 } terms on the array diagonal (row k 
and column k) and each E{S£s£ *} term at a respective row position k and column position 
m, this set of values is known as the "spatial covariance matrix". The entire process may be 
carried out independently for each DFT frequency.bin, i.e. each value of transformer output 
5 index n. 

Using some or all of the covariance terms in the spatial covariance matrix, the bearing (or 
angle of incidence or arrival) of one or more received signals may be estimated using 
standard techniques: see e.g. H L Van Trees "Optimum Array Processing" (part IV of 
"Detection, Estimation and Modulation Theory"), Wiley, New York, 2002, which discloses for 
10 example the MUSIC algorithm and least squares fitting. This is well known in the art of 
direction finding and will not be described. 

Signal power received by antennas 12 may fluctuate: because of this it is advantageous to 
carry out averaging in connection with the measured quantities P^, and P^b defined 
above. It is convenient to define a "commutation cycle" as consisting of a cycle of collecting 

15 data as described above from each of the settings b to e, p to r, of the three switches SW1 , 
SW2 and SW3 required to determine P^, P^ and P^b for one value of the frequency bin 
index n. Further data is then collected using additional commutation cycles, each cycle 
giving a respective set of values of P,^ and P^b. The resulting multiple values of P kk , 
PkmA and Pkms are then used to provide average values of each. These averages are then 

20 used in Equations (1) to (1 0) above. 

The rate at which commutation cycles are carried out is referred to herein as the 
"commutation cycle rate". One or more of the signals received by the antennas 12 may 
contain periodic fluctuations, if for example data symbols are carried by a signal at a 
particular rate. If the period of the fluctuations is equal or close to a multiple or sub-multiple 
25 of the commutation cycle rate, then covariance terms derived as described above may be 
subject to consistent or systematic errors (known as "biases"). These may result in errors in 
estimated angles of incidence. To avoid this problem, the order in which the commutation 
(switch setting) is performed within each commutation cycle may advantageously be varied 
in a random or pseudo-random sequence between successive cycles. 
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Equations (1) to (10) above may be re-evaluated for other values of the frequency bin index 
n, e.g. adjacent frequency bands corresponding to index n-1 and n+1. Output bearing 
estimates are then obtained for various values of n and are combined to give an estimate 
for which error limits can be calculated. 

5 A simulation was made of a direction finding system 1 0 of the invention having four 

antennas 12 in a square array. The simulation .envisaged two emitters with identical carrier 
frequencies to be located in bearing in the plane of the antenna array. The emitters were 
treated as both transmitting at 6.25ksymbol/sec with a modulation type of QPSK, as 
described by R E Ziemer and R L Peterson, "Introduction to Digital Communication", 
10 Maxwell Macmillan International, New York, 1 992. They were at respective bearings' or 
angles of arrival of 30° and 70° at the antenna array relative to a predefined reference 
direction. The signal to noise ratio (emitted signal power divided by total noise power over 
25.6MHz bandwidth) was assumed to be high. The sampling rate was 25.6MHz and a 
block of 4096 samples were collected, a Hamming window was applied (see Oppenheim et 
al. mentioned above), and the result was processed by discrete Fourier transform. The 
sample collection time was 0.16 milliseconds, which corresponds to one emitter symbol 
period. If only one block of samples were to be processed, results are poor, because 
signals are highly correlated over the sample collection time. First and second blocks of 
samples were collected with an intervening time interval corresponding to at least three 
symbols to avoid this. Emitters of interest often contain a root raised cosine (RRC) pulse 
shaping filter which is shorter than this intervening time interval, in which case there is no 
correlation between signals in the sample blocks. The overall covariance matrix was then 
derived using parameters obtained' using Equations (1) to (10) for both data blocks and 
averaging, Matrix elements were processed using the known MUSIC algorithm (see Van 
Trees mentioned above) to determine emitter bearings. This algorithm is neither the only 
nor even the best algorithm for this purpose, but it is simple and often suggested for that 
reason. 

Figure 2 is a graph of MUSIC spectrum against angle of arrival (AoA) derived from 
covariance matrix elements obtained using the method of the invention. It shows 
30 reasonably well-defined peaks 40 and 42 at angles of arrival of 30° and 70°, showing that 
the simulated emitters have been located. 
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For comparison, a trial was made of the single receiver technique disclosed byCMS See, 
"High Resolution DF with a Single Channel Receiver", 0-7803-7011-2/01 IEEE 2001. Using 
the same parameters and MUSIC algorithm processing, this technique gave the results 
shown in Figure 3. Here the emitter at 30 degrees has not been resolved, arid the 70 
5 degree emitter has been resolved at 50, but with a peak height ~15x smaller than the peak 
42 (NB Figures 2 and 3 have axes with different scales).- 

For completeness the conventional method (which requires one receiver per antenna) was 
also simulated using the same parameters and MUSIC algorithm processing. This method 
gave the results shown in Figure 4. It shows very well defined peaks 60 and 62 at angles of 
10 arrival of 30° and 70°, showing that the simulated emitters have been located with better 
accuracy than the invention, but this is at the expense of using one receiver per antenna 
instead of one receiver only for all antennas. 
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Claims 

1 . A direction finding system incorporating a plurality of antennas (12) characterised in 
that the system (10) also includes: 

a) means (SW1) for determining individual antenna signal strengths; 

b) combining means (SW2, P2, P3, SW3, 18) for deriving combined antenna 
signal strengths by forming combinations of first and second antenna signals 
derived from different antennas (12), wherein the second antenna signals are 
in two sets with signals in one set having a non-zero phase difference relative 
to signals the other set; and 

c) means for deriving covariance matrix elements from antenna signal strengths 
and for determining at least one emitter bearing therefrom. 

2. A direction finding system according to Claim 1 characterised in that the relative 
phase difference is in the range 30 to 120 degrees, and the means (SW1) for 
determining individual antenna signal strengths and the combining means (SW2, P2, 
P3, SW3, 18) are arranged to enable successive signal strengths to be derived in 
successive steps. 

3. A direction finding system according to Claim 2 characterised in that the relative 
phase difference is substantially 90 degrees. 

4. A direction finding system according to Claim 3 characterised in that the combining 
means (SW2, P2, P3, SW3, 1 8) is arranged to combine antenna signals with equal 
gain magnitude and with equal and unequal phase. 

5. A direction finding system according to Claim 1 characterised in that the combining 
means incorporates phase shifting means (P3) switchable into and out of an 
antenna signal path. 
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6. A direction finding system according to Claim 1 characterised in that the combining 
means incorporates an adder (1 8) having two inputs both switchably connected to 
individual signal paths extending to respective antennas (12). 

7. A direction finding system according to Claim 1 characterised in that: 

a) the means for determining individual antenna signal strengths comprises a 
first multipole switch (SW1) having input poles (b, c, d, e) connected to 
receive signals from respective antennas (12); 

b) the combining means (SW2, P2, P3, SW3, 18) incorporates a second 
multipole switch (SW2) having input poles (b, c, d, e) connected to receive 
signals from respective antennas (12) and a third multipole switch (SW3) for 
switching phase shifting means (P3) into and out of an antenna signal path 
extending via the second multipole switch (SW2); and 

c) the combining means (SW2, P2, P3, SW3, 18) also incorporates adding 
means (18) for combining signals, the adding means being arranged to add 
an antenna signal in a first signal path extending via the first multipole switch 
(SW1) to another antenna signal in a second signal path extending via the 
second and third multipole switches (SW2, SW3). 

8. A method of direction finding using a plurality of antennas (12) characterised in that 
the method incorporates: 

a) determining individual antenna signal strengths; 

b) determining combined antenna signal strengths by forming combinations of 
first and second antenna signals derived from different antennas (12), 
wherein the second antenna signals are in two sets with signals in one set 
having a non-zero phase difference relative to signals the other set; and 

c) deriving covariance matrix elements from antenna -signal strengths and 
determining at least one emitter bearing therefrom. 
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A method according to Claim 8 characterised in that the relative phase difference is 
in the range 30 to 120 degrees and successive signal strengths are determined in 
successive steps. 

A method according to Claim 9 characterised in that the relative phase difference is 
substantially 90 degrees. 

A method according to Claim 10 characterised in that the step of forming combined 
antenna signal strengths combines antenna signals with equal gain magnitude and 
with equal and unequal phase. 

A method according to Claim 8 characterised in that the step of forming combined 
antenna signal strengths includes switching phase shifting means (P3) into and out 
of an antenna signal path. 

A method according to Claim 8 characterised in that the step of forming combined 
antenna signal strengths includes adding signals in signal paths extending 
switchably to different antennas (1 2). 

14. A method according to Claim 8 characterised in that: 

a) the step of determining individual antenna signal strengths comprises 
switching signals from antennas (12) via a first path; 

b) the step of forming combined antenna signal strengths incorporates: 

0 switching signals from antennas (1 2) via a first path for combining; 

ii) switching signals from antennas (12) via a switch selectable second 
path or a third path for combining, the third path being arranged to 
phase shift antenna signals therein relative to antenna signals in the' 
second path; and 

iii) adding a first path antenna signal to second and third path antenna 
signals individually. 
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ABSTRACT 



Direction finding by radio comprises arranging an array of antennas (12) to receive signals 
from emitters, selecting individual antenna signals using a first multipole switch (SW1) and 

5 determining antenna signal strengths. Individual antenna signals are also selected by a 
second multipole switch (SW2), which routes a selected signal to a third multipole switch 
(SW3). The third switch (SW3) switches a phase shifter (P3) into and out of an antenna 
signal path. An adder (18) is employed to add an antenna signal in a first signal path 
extending via the first multipole switch (SW1) to a different antenna signal in a second 

10 signal path extending via the second and third switches (SW2, SW3). This determines 
combined signal strengths between pairs of antenna signals, one of which either has or has 
not been relatively phase shifted depending on the third switch position. Covariance matrix 
elements are determined from signal strengths enabling emitter bearings to be derived. 
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